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In this paper we show that the claims in Jjj related to our analysis in Q are wrong. 



As is well known, current observational evidence strongly favors an accelerating and spatially flat Friedmann- 
Robertson- Walker universe || , Q . Since normal matter obeys the strong energy condition and cannot drive acceler- 
ated expansion, recourse is often made either to a small cosmological constant or an almost evenly distributed source 
of energy -termed "quintessence"- in the form of a self-interacting scalar field with equation of state p^ = (7^ — l)/?^ 
(where < 7^ < 1) such that it provides a negative pressure high enough to render the deceleration parameter 
negative (i.e., q = —a/(aH 2 ) < 0), see e.g. ||. As usual, a denotes the scale factor of the FRW metric and H = a/a 
the Hubble factor. Owing to the fact that in general the scalar field and matter energy are expected to scale dif- 
ferently with expansion, the question arises: "why the ratio between matter and quintessence energies should be of 
the same order precisely today?" Put another way, "where the relationship (^i m /^<t>)o — O(l) comes from?" This 
is the so-called coincidence problem [^J. Here the zero subindex means present time and the dimensionless density 
parameters of matter and quintessence arc defined by 

n» - - (« = = i), 

respectively. Note that the cosmological constant case is recovered by setting 7^ = 0. 

Obviously, solving the coincidence problem in an accelerating universe amounts to show that the ratio Q, m /Q,^, tends 
to a constant for large times with q < 0. Then, one may choose the free parameters of the model to fix the constant 
ratio to order unity. In a recent paper we demonstrated that the coincidence problem and an accelerated expansion 
phase of FRW cosmologies cannot be simultaneously addressed simply by the combined effect of a perfect matter fluid 
and a quintessence scalar field. Nonetheless, if the matter fluid is not perfect but dissipative, both problems can find 
a simultaneous solution for FRW spatially flat and open FRW universes . 

A key point in our derivation was the stationary condition 

TT 2H 

7m + — =70 = -XT72, (1) 
Pm O-TZ 

as well as the violation of the strong energy condition 

TT < (| ~ 7m) Pra (2) 

(cf. ||), where, as usual, 7 m denotes the baryotropic index of matter, i.e., p m — (j m — l)p m thereby 1 < j m < 2. 
Likewise it stands for the dissipative stress which is negative for expanding universes. It may be associated with 
particle production JtJ, may be understood as frictional effects arising in mixtures Q or even model other kinds of 
sources (e.g. a string dominated universe as described by Turok or a scalar field). 

Equation (|l|) expresses the condition for fi m and Q<p to be constants. It can be straightforwardly derived by com- 
bining the conservation equations for the matter fluid and quintessence field 
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Pm + 3 ( 7 m + — ) p m H = 0, and p^ + Sj^p^H = 0, (3) 

V Pm) 

with the definitions of the density parameters 17 m and f2^, respectively. Thus, the condition fl m = Cl^ = translates 
into @. Equation (|) expresses the asymptotic stability of the solution Q = Q m + f2^ = 1. By slightly perturbing 
the solutions f2j — constant (i = m, <fi) one finds they are stable for spatially flat and open accelerating universes -see 
Eqs (18) and (22) of §]. At this stage we would like to emphasize two points, namely: (i) our derivation shows that a 
dissipative stress negative enough to satisfy the Eqs. (|l]) and (^) is required to provide a coincidence solving attractor 
solution. This derivation does not hinge on the specific choice of the potential of quintessence field, (ii) Our proof is 
meant just for universes under accelerated expansion, i.e., for q = |7^ — 1 < 0. Clearly, the conditions expressed by 
Eqs. (|l) and (||) define the set of transport equations for 7r or, equivalently, the set of expressions of the bulk viscosity 
coefficient, leading to accelerated stable solutions (see §IV of for some examples). 

However, recently Ibahez et al. studied the autonomous system of equations of a FRW universe filled with a dissipa- 
tive matter fluid and a self-interacting scalar field V(4>) oc exp(k(j>) and found stable and unstable equilibrium points 
with f2j = constant [Q . The stability of these points seems to depend on the specific equation of state of ir as well as 
on the values assumed by different parameters of their model. Since this outcome looks at variance with point (i) of 
the precedent paragraph Ibahez et al. erroneously claimed to have found examples that contradict our findings. The 
fact is, however, they overlooked point (ii): all the unstable solutions in p| correspond to non-accelerating universes. 
To be specific: the equilibrium point of Eqs. (29) has q > (cfr. Eq. (32)). The set of equilibrium points associated 
to a massless scalar field (r = 0) have also q > 0. Likewise, the equilibrium point of Eqs. (40)-(42) is either stable and 
accelerating or unstable and decelerating (see Eq. (48) and Fig. 4). Again, the equilibrium point of Eqs. (49)-(52) is 
non-accelerating (see Eq. (53)). Finally, the equilibrium point (54)-(56) has q > 0. 

Therefore they do not invalidate in any way whatsoever the findings of 0] . We believe this is more than enough to 
dismiss the claims of Ibahez et al. However, we have found inconsistencies in their analysis. In the remaining of this 
paper we concentrate in bringing them to the fore. 

The starting equations of Ibahez et al. are: 

B=-B 2 - X - (p m + 3p m + 3tt + 20 2 - 2V{4>)) (4) 

1 • "W( 

3H 2 = Pm + -4? + V(4>) ~— {K= ±1, 0) (5) 

2 a z 

p rn = -3H(p m + p m + 7f) (6) 

i-^Hi-Z&. (7) 

The dissipative pressure tt is assumed to satisfy the truncated Israel-Stewart equation |To| 

n + Tir = -3Cff, (8) 

where £ is the coefficient of bulk viscosity and r the relaxation time (C > 0, r > 0). Strictly speaking, Eq. Q is 
valid only when the fluid is close to equilibrium. However, let us assume it holds even when the fluid is far from 
equilibrium. (A more rigorous and comprehensive analysis would make use of the full transport equation of the 
Israel- Stewart theory.) Ibahez et al. also assumed the linear baryotropic equation of state p m — ("f m — l)p m , but with 
7 m = constant, and two different relations for the coefficient of bulk viscosity and relaxation times. To translate 
(3)-(|8|) into an autonomous system of differential equations Ibahez et al. introduced the set of variables 

3H 2 ' H 2 V6i? 

3 H 2 ' T 2' V ; 



as well as a new time parameter r defined by 



dr = H(t) dt. (10) 
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The Friedmann equation (4) now reads 



i - ni - * 2 



K 



(11) 



and the autonomous system takes the form 
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(12) 
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(14) 

(15) 
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where the prime denotes derivative with respect to t, and x = 1 + ^(37 m — 2)f2 m + + 2\I> 2 — T. 

Note that Eq. ( |i~3| ) is the truncated Israel-Stewart transport equation, while Eq. p6| ) is just the derivative of the 
change of variables introduced in @ . In Eq. ([b|) , first case of Ibanez ef a/. , use has been made of the relationships 
C = ocpm> T = a Prn l i introduced by Belinskii et al. plj . In their second case a somewhat modified relationships were 
used for £ and r (introduced by their Eqs.(37)), amounting to replace ( |l3[ ) by a similar equation. 

The equilibrium points are found by setting Q! m = £' = iff' = V = h' = 0. As a consequence, p m , 7T, 0, V^($>) 
and H are constants there (use of the variables introduced in Eqs. (||) and (|l(]) excludes H — 0). Then, in virtue of 
Eq. (0) we have dV (<fi) / d(f> — and = 0, i.e., V(^) must have an extremum at the equilibrium points, and 7^ = 0. 
Potentials having this feature yield de Sitter solutions on the equilibrium points and for them, as noted above, we 
could replace the quintessential field with an effective cosmological constant. 

From Eq. (4), with H, p, V(4>) constants and (f> = 0, we get K = 0. Moreover, from (5) it follows that ir = —j m Pm- 
This is independent of the specific form of V{<j>). Notice that this is a particular case of the result already found in 
0, as follows from (|l[) and the assumption H = constant. 

The setting h! = by Ibahez et al. (their Eq.(18)) leads to wrong conclusions. In the first case (Belinskii et al. 
relationships), they erroneously state that its solution is h = rather than h = — 3a(3fi m ) n /£ (with x — 0), as follows 
from their Eqs. (15) and (19); and h — constant in the second case. Besides, their relationship h' = neither belongs 
to the set of Einstein equations nor describes any property of the sources of the gravitational field. 

Altogether, the work of Ibanez et al. misinterpret our results (to say the least) rendering the claims in Q void. 
Further, even if the analysis of Ibahez et al. were correct (that is not), it would fail to disqualify our findings in Q as 
all their counterexamples correspond to non-accelerated cosmic expansions. 
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